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In a class of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, we derive both vacuum and inte¬ 
rior Schwarzschild solutions under the condition that the derivatives of a scalar field (j> with respect 
to the radius r vanish. If the parameter oh characterizing the deviation from Horndeski theories 
approaches a non-zero constant at the center of a spherically symmetric body, we find that the 
conical singularity arises at r = 0 with the Ricci scalar given by R = —2a.Yilr^. This originates 
from violation of the geometrical structure of four-dimensional curvature quantities. The conical 
singularity can disappear for the models in which the parameter an vanishes in the limit that r —>■ 0. 
We propose explicit models without the conical singularity by properly designing the classical La- 
grangian in such a way that the main contribution to an comes from the field derivative (/>'(r) 
around r = 0. We show that the extension of covariant Galileons with a diatonic coupling allows for 
the recovery of general relativistic behavior inside a so-called Vainshtein radius. In this case, both 
the propagation of a fifth force and the deviation from Horndeski theories are suppressed outside a 
compact body in such a way that the model is compatible with local gravity experiments inside the 
solar system. 


I. INTRODUCTION 

Over the past few decades, there have been numerous 
attempts for extending General Relativity (GR) to more 
general gravitational theories [I|. This is motivated by 
the ultra-violet completion of gravity ii or by the ob¬ 
servational evidence of early and late phases of cosmic 
acceleration [ 1 , In particular, the problem of dark en¬ 
ergy implies that the today’s acceleration of the Universe 
may be related to some infra-red modification of gravity. 

If we extend GR in such a way that only one addi¬ 
tional scalar degree of freedom (DOF) arises, Horndeski 
theories @ are known as the most general scalar-tensor 
theories with second-order equations of motion Q. The¬ 
ories with derivatives higher than second order can be 
prone to instabilities due to the appearance of extra ra¬ 
diative DOF 0. Nevertheless, it is possible to generalize 
Horndeski theories to those with only one propagating 
scalar DOF. In Hofava-Lifshitz gravity Q, for example, 
addition of higher-order spatial derivatives does not in¬ 
crease the number of the scalar DOF. 

Another possibility for the extension of second-order 
gravitational theories is to violate two conditions Horn¬ 
deski theories obey One of such conditions is 

related to a geometric modification of the Einstein- 
Hilbert Lagrangian Leh = M^^R/2, where R is the 
four-dimensional Ricci scalar and Mpi is the reduced 
Planck mass. In terms of the 3-I-I Arnowitt-Deser-Misner 
(ADM) decomposition of space-time 0 , this Lagrangian 
can be expressed in the form Leh = Ai {K'^ — Kf„yK^'^)-\- 
BifR, with —A 4 = B 4 = Mpj/2, where K and TZ are the 
traces of three-dimensional extrinsic and intrinsic curva¬ 
tures respectively. On the isotropic cosmological back¬ 
ground Horndeski theories contain the same form of La¬ 
grangian as Leh, but A 4 and B 4 are functions of the 


scalar field (f and its kinetic term X. These functions 
satisfy the particular relation A 4 = 2 XB 4 X — B 4 [llj, 
where B 4 X = dB 4 ldX. 

GLPV theories do not obey this particular relation as 
well as another condition associated with the Horndeski 
Lagrangian L 5 Q. In the unitary gauge, where the con¬ 
stant field hypersurfaces are identified with the constant 
time hypersurfaces, the Hamiltonian analysis shows that 
GLPV theories possess only one scalar DOF [l^. The 
simplest class of GLPV theories corresponds to functions 
A 4 and B 4 depending on <f) alone, in which case the dif¬ 
ference between —A 4 and B 4 characterizes the deviation 
from Horndeski theories. This deviation gives rise to sev¬ 
eral interesting observational signatures such as the mix¬ 
ing of scalar and matter propagation speeds 
the large gravitational slip parameter jl5l| , and the real¬ 
ization of weak gravity |16l| . 

For the consistency with local gravity experiments, 
the interaction with matter mediated by the scalar DOF 
should be suppressed inside the solar system m. In 
Horndeski theories, nonlinear scalar-field self interactions 
can screen the fifth force [3 M through the Vain¬ 
shtein mechanism [ 2 ^. In GLPV theories it was recently 
claimed that the Vainshtein mechanism tends to break 
down around the vicinity of a compact object [3 - [3| . 
Whether this conclusion holds or not for all classes of 
GLPV theories remains to be an open question. 

In this paper, we revisit the analysis of spherically sym¬ 
metric solutions in GLPV theories to study how much de¬ 
viation from Horndeski theories can be allowed. On using 
the background equations of motion derived in Sec. |lTl 
we first obtain vacuum solutions around a point source 
(r = 0) in Sec. IHII and show that the Ricci scalar is 
given by i? = — 2aH/r^ in the absence of the cosmolog¬ 
ical constant A. Since an 7 ^ 0 in GLPV theories, there 
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is a conical singularity at the origin. In Sec. IIVI we also 
derive solutions inside a compact object for constant an 
under the condition that radial derivatives of van¬ 
ish everywhere and show that R is again divergent at 
r = 0. Since the boundary condition = 0 at r = 0 
is most natural for the regularity of solutions, the conical 
singularity inevitably arises for the theories in which an 
approaches a non-zero constant as r —> 0. 

A possible way out for avoiding the conical singularity 
is to consider the theories with a vanishing an at the 
origin by designing the Lagrangian such that the main 
contribution to an corresponds to A-dependent terms. 
The exact form of choice for the Lagrangian is essential 
to remove such a conical singularity, and the absence of a 
symmetry which preserves its form might be a problem, 
once one considers, e.g., quantum corrections. In Sec. m 
we discuss conditions for eliminating the conical singu¬ 
larity in the framework of GLPV theories and present 
an explicit Lagrangian which allows for this possibility. 
In Sec. EH we derive the field profile and gravitational 
potentials inside and outside the body under the approx¬ 
imation of weak gravity for a concrete model without the 
conical singularity. We show that the Vainshtein mech¬ 
anism is efficient enough to suppress the fifth force and 
the deviation from Horndeski theories. We conclude in 

Sec. Em 


II. BACKGROUND EQUATIONS OF MOTION 


We consider the theories described by the action @ 

4 

s= f f , ( 2 . 1 ) 


where g is the determinant of metric and is 
the Lagrangian of matter fields 'i’m with the energy- 
momentum tensor Tjf = diag{—pm, Pm, Pm, Pm)- We 
assume that matter is minimally coupled to the metric 
The Lagrangians L 2 ,L^, and L4 are given by 


L 2 

L 3 

Li 


A2(</),A), (2.2) 

(C'3 + 2AC'3,x)n</> + AC'3.,^, (2.3) 

BiR - [(□</,)2 _ 


2 {Bi + A 4 - 2XBi^x) 




A 2 -r- -r ^ 

, (2.4) 


where denotes covariant derivatives, A 2 ,C 3 , A 4 , B 4 
are functions of the scalar field (j) and its kinetic energy 
X = g^'^X Introducing the function F4 defined 
by 


- FiX^ = B 4 


A4 - 2 XB. 


4,X , 


(2.5) 


one can write the Lagrangian L 4 of the form 211 


L4 = B4R - 2 B 4 ,x [(□</>)" - 


where t^upa is the totally antisymmetric Levi-Civita ten¬ 
sor. Horndeski theories satisfy the condition F 4 = 0, i.e., 
B4 + A4 — 2XB4 ^x = 0, under which the terms after the 
second lines of Eq. (1^ vanish. 

The full action of GLPV theories contains the La¬ 
grangian L5 associated with the Einstein tensor, but we 
do not take into account such a contribution in this pa¬ 
per. In fact, inclusion of L5 tends to prevent the recov¬ 
ery of GR in the solar system even in Horndeski theories 

Using the 3-1-1 ADM formalism, we can construct sev¬ 
eral scalar quantities such as K = and TZ = 

g^^TZpu from the extrinsic curvature and the intrin¬ 
sic curvature TZf^i, [ll|. In the unitary gauge {(p = (pit)), 
the Lagrangian L = L 2 F L 3 + L 4 is equivalent to 
L = A 2 +MK + AiiK"^ - -t Bill with the 

relation 


A3=2|A|3/2(^C'3,x + ^) , (2.7) 


where the sign of X is different depending on the given 
space-time. On the spherically symmetric and static 
background we have that A > 0, whereas on the isotropic 
and homogenous background A < 0. 

The deviation from Horndeski theories can be quanti¬ 
fied by the parameter Q 


2XB4 X-B4 FiX^ 
an = -^-1 = 


A 4 


A 4 


( 2 . 8 ) 


We also define 


at = - 


A 4 


- 1 . 


(2.9) 


For linear perturbations on the spherically symmetric 
background, at is related to the tensor propagation speed 
squared c^, as at = l/c? — 1 [H, [l^. Then, at charac¬ 
terizes the deviation of Ct from 1. If the function B4 
depends on (p alone, an is equivalent to at. 

We consider the spherically symmetric background de¬ 
scribed by the metric 


ds^ = -t r^{de^ + sia^ 9 dip^) , 

( 2 . 10 ) 

where the gravitational potentials 4'(r) and $(r) are 
functions of the distance r from the center of symme¬ 
try. Variation of the action m leads to the equations 
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of motion 

'4e-2‘fA4 


— Cl + 4C2 ) — A2 + C3 — C4 


2A4 , _ 24 > 1 

(e -l-at)=-p, 
4e-2'fA 


— Cl + 4C2 I 4/^ + A2 — 2A^A2 X — 


+ — [^ 4(0 — 1 — q;h) + »’C 2 ] — —Pn 


( 2 . 11 ) 

2Ci 

r 

( 2 . 12 ) 


2e-2-J>A4 ( 4 '" + 


2 e-‘*’A 4 1 . , 

—— + -C 4 r ) 4.' 


Cl - 2 


6 - 2 -^ A 4 


+ C 2 (l + r4'') 




+A2 — C3 + —C4 — —Pm I 
Pm + {Pm + Pm) = 0 j 


(2.13) 

(2.14) 


where X = e a prime represents a derivative with 

respect to r, and 


Cl = 2e-^XA3.3c , C 2 = 

r 

C 3 = (A3_0 + 2e“^^0"A3_x) , 

^ {A4,4, + 2e-^^rA4,x) 


4,X 


(2.15) 


In GR we have —A 4 = B 4 = M'^J2 and an = at = 
0. In GLPV theories, the parameter an gives rise to a 
contribution comparable to the dominant term A 44 >/r^ 
in Eq. (12.121) for |an | larger than the order of 14)|. 

In the following we will assume that both A 4 and B 4 
also remain finite as r —0. For this class of theories 
(without shift symmetry), on a static spherically sym¬ 
metric background, the field (j) is required to be only 
r-dependent, and 4 >'(r ^ 0 ) =0 for the regularity in 
compact objects. In this case, we have that X —0 at 
the origin. 

It should be noticed that the action (EH) remains finite 
in the limit X —?► 0. In fact, for analytic profiles for 
the field, one can verify that 1/—^^4 —?► constant around 
the origin. This should not come as a surprise as, in 
fact, the same limit is well-defined also for Horndeski 
theories, whose Lagrangian reduces to the first line of 
Eq. (I2.4E We also note that the equations of motion 
(I2.11I) - (I2.13I) do not contain the X-dependent term in 
the denominators. It is clear that finite values of A 4 and 
B 4 are allowed as long as $',4'' —0, as requested for 
standard boundary conditions for compact objects. 


III. VACUUM SCHWARZSCHILD SOLUTIONS 

First of all, we derive exact solutions to Eqs. (I2.11I) 
and (I2.12P in the absence of matter for a point source 
with mass Xi. Since the field derivatives ()>'(r) and ()>"(r) 


vanish in this case, the terms Ci defined in Eq. (I2.15I) do 
not contribute to the equations of motion. The terms 
A 4 and B4 are regarded as constants, with the particular 
relation an = ott = —B4/A4 — 1 . The term A 2 is related 
to the cosmological constant A (> 0), as A 2 = —A. Then, 
Eqs. (12.111) and (12.121) reduce to 


4e-^‘^A4 


4>' + A-^(e- 2 ®-l-aH) = 0 j( 3 . 1 ) 


4 e 


r 

- 2 ^ A 


4 /' - A + ^(e-2* - 1 - an) = 0 .(3.2) 


r 

The solutions to these equations are given by 


= ( 1 -I- an 


24/ 


A 


Cl 

2 A 4 r 
A 


= - 6 A 4 C 2 1 1 -I- an + -^1" + 


Cl 

2 A 4 r 


(3.3) 

(3.4) 


where Ci and C 2 are integration constants. The devia¬ 
tion from Horndeski theories leads to the contribution 
to gravitational potentials. The cosmological constant 
gives rise to the contribution Kr^/{QA 4 ), as it happens 
for the Schwarzschild-de Sitter solution. The two coef¬ 
ficients Cl and C 2 also arise in GR. The constant ci is 
chosen as ci = — 4 A 4 GA 4 to obtain the standard term 
—2GA4/r in Eqs. (j3.3|) and ()3.4I) . To recover the usual 
relation the constant C 2 is chosen to be 

C 2 = — 17 ( 6 ^ 4 ). In this way, we have fixed two freedoms 
(Schwarzschild mass Ai and the time reparametrization) 
as in the context of GR. 

Then, we obtain the following solution 


^24- ^ 


1 -f an + 


A 


2GM 

r 


(3.5) 


This is the extension of the Schwarzschild-de Sit¬ 
ter solution with the additional factor an- Since 
g**Rtt = g^'^Rrr = -A/( 2 A 4 ) and g^^Res = g^^R^^ = 
—K/{ 2 A 4 ) — a-a/r^, the Ricci scalar is given by 


R=- 


2A 


2 aH 


(3.6) 


which diverges at r = 0 for an 7 ^ 0. Provided that 
an 7 ^ 0 for r —>■ 0, the divergence of R occurs independent 
of the choices of ci and C 2 , so it cannot be eliminated by 
the change of boundary conditions. 

The divergence of curvature at r = 0 can be in¬ 
terpreted as a conical singularity originating from the 
0, ip contributions to R. This singularity cannot be 
eliminated even for A —7 0 and A4 —7 0. In this 
limit, the three-dimensional spatial metric is given by 
^*( 3 ) = (1 + ctR)~^dr'^ + r^(c?0^ -f sin^0d(p^). Defin¬ 
ing f = r/-\- an, the two-dimensional metric in the 
0 = 7r/2 plane is of the form = c?f^ -|- , where 

(p = 1/1 + o^H Since the angle ip is not restricted to 
be between 0 and 27r, the resulting space-time is a cone 
with a geodesically incomplete point at r = 0. 
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To avoid such a conical singularity in the vacuum, we 
need to impose the condition ajj =0. In the next section 
we will show that the same conclusion also holds in the 
presence of a matter source for constant an- 


IV. INTERIOR SCHWARZSCHILD SOLUTIONS 
FOR CONSTANT an 


We derive interior and exterior Schwarzschild solutions 
under the conditions that the density Pm inside a com¬ 
pact body is constant and that (j)'{r) = 0 and = 0 

everywhere. Again, we deal with an as a constant satis¬ 
fying the relation an = at = —-B 4 /A 4 — 1 . 

Integration of Eq. (12.141) gives 

Prn = -Pm+Be-^, (4.1) 


where S is a constant. Solving Eq. (I2.11|) for $, we obtain 
the same solution as Eq. (13.31) with the replacement A —>■ 
A + Pm ■ For the regularity of metric at r = 0 we need to 
choose Cl = 0 , so the solution reduces to 


= (1 -I- an - ^ , 


where 


_ A -|- Pm 

6A4 


(4.2) 


(4.3) 


There exists the term an in Eq. (03, which cannot be 
eliminated by any boundary condition. As we will see 
below, this again leads to the conical singularity at r = 0. 
Solving Eq. (12.121) for T, we obtain 


——- --V^l + an-Ar'^ , (4.4) 

2(A -|- Pm) 


where V is a constant. There are two free integration con¬ 
stants B and V in addition to the mass A4 of the compact 
object. These constants can be fixed by matching inte¬ 
rior and exterior solutions of $ and at the radius rg of 
the body determined by the condition Pm{fg) = 0. Anal¬ 
ogous to the derivation of Eq. (1531) . the exterior vacuum 
solution for r > rp is given by 

24/ —-I ^ A ‘2 2GA1 / i r' \ 

= e = 1 +an - - Ar^ -. (4.5) 

A + Pm r 

Matching the solutions at r = tq, it follows that 


B = Pmy'T+'aH— 

^ _ 2GA4 A + Pm 

.A o 5 

^0 Pm 

jy Pm — 2A 

2(A -|- Pm) 


(4.6) 

(4.7) 

(4.8) 


After substitutions of Eqs. (14.61) - (14.81) into Eqs. (14.21) 
and (HU), we obtain the interior Schwarzschild solution 


with appropriate matching conditions. Then, we can find 
the behavior of three curvature scalars R, S = , 

and T = R^vXpR^'^^'^■ In the limit where r —>■ 0, it follows 
that 


R 


o rp 4a^ 

^2 ’ ^4 ’ ^ J .4 ■ 


(4.9) 


We have dropped next-to-leading order corrections, 
which appear as a constant for R and as the terms in¬ 
volving an!A for S and T. Taking into account these 
corrections, the Gauss-Bonnet term Rq^ = R"^ — AS + T 
is proportional to an/‘A around r = 0. The integration 
constants have been fixed as Eqs. (14.61) - (I4.8p . but the 
singular behavior of i?. S', T is independent of the choice 
of boundary conditions. Thus, for an = constant, the 
conical singularity similar to that discussed in Sec. Illll is 
present at the center of a compact object. 

It should be noted that the metric and matter fields 
are all regular at the origin, whereas the singularity ap¬ 
pears only at the level of the Ricci scalar, i.e., for tidal 
forces and geodesics deviation. Although we have found 
the presence of this singularity for exact solutions (which 
always give a clear insight of theories), one may worry 
that it is merely a consequence of the fact we have as¬ 
sumed a specific constant field profile (j)(r) = constant 
everywhere. But this “special” solution is actually “the” 
solution which is approached at the origin for compact 
objects (assuming the analytic profile for each field). 

In fact, if the fields are analytic around r = 0, then we 
can Taylor-expand them as: (j) = (j)c + 12 -|- 0{A), 

$ = 4>c + C 4 > ?'^/2 + 0 (?'^), 4/ = Aic + A/ 2 -\- 0 {A), and 
p = Pc + Cp A12 + 0{r^). The coefficients ('s are numer¬ 
ical constants which can be determined by the equations 
of motion or boundary conditions. Here, the crucial point 
to be shown is whether the variation of (j){r) induced by 
the term C,^A /2 can modify the existence of the conical 
singularity. 

We do not specify the functional forms of A 2 , ^ 3 ,^ 4 , 
but the leading-order contributions to the terms A 2 , 
As^x, A 3 _ 0 , A 4 ^x and are regarded as constants. Us¬ 
ing the above Taylor-expansions of (j) and $, the leading- 
order terms in Eq. (I2.15|) are given, respectively, by 
Cl = cir^, C 2 = C 2 r, C 3 = csr, and C 4 = £ 4 , where cCs are 
constants associated with In the limit that r —>■ 0 the 
terms Ci, C 2 , and C 3 vanish, so only the non-vanishing 
contributions to Eqs. (I2.11D - (I2.13I) are —C 4 in Eq. (12.111) 
and C 4/2 in Eq. (12.131) . These C 4 terms behave in the 
same way as another constant A 2 = —A already studied 
above. This discussion shows that the variation of (j){r) 
around r = 0 does not give rise to any new contribution 
to modify the existence of the conical singularity. 

More concretely, one can verify from the integration 
of Eq. (12.111) that the gravitational potential e~^'^ ac¬ 
quires the terms proportional to cir®, C 2 r^, cgr^, C 4 A, 
and in addition to those given in Eq. (14.5|) . In the 
limit that r —>■ 0 all these additional terms vanish, so 
the gravitational potential again approaches the value 
—>■ 1 -(-aH- Similarly the additional terms Ci and C 2 
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appearing in Eq. (12.121) only give rise to vanishing con¬ 
tributions to 'I' as r —>• 0, such that —>■ 1 -I- an for 

r —>■ 0. In addition to the fact that the integrated solution 
631) of Eq. (12.141) is used in this procedure, the above 
solutions are consistent with Eq. (12.131) . Hence they are 
the solutions of the full equations of motion (I2.11I) - (I2.14I) 
around r = 0. The existence of the term an in $ and d' 
shows that the variation of the field (j){r) around r = 0 
cannot eliminate the conical singularity. 

The robustness of the behavior = 1 -I- an 

for r —)► 0 can be also confirmed by directly substituting 
the above expansions into the background equations of 
motion. On expanding Eq. (12.111) around the origin, we 
find that the first non-trivial contribution which needs to 
be canceled is given by 

4 0) + e-2^'= A4{cj)c, 0)], (4.10) 

from which we obtain 6“^^° =1-1- q;h(<(>c,0). Applying 
the same procedure to Eq. (12.121) with Eq. (14.11) , it follows 
that 6 ^'^'= = l-|-aH(</'cj 0). Thus, the theories approaching 
a non-zero value of an as r —> 0 lead to the appearance 
of the conical singularity. 

It should be noticed that we have focused on static so¬ 
lutions and that we have not imposed the shift symmetry 
of the theory. Abandoning any of the two assumptions 
could lead to a different behavior for the theory, and pos¬ 
sibly, to the disappearance of the singularity. We leave 
this interesting topic for a future study. 

V. CONDITIONS FOR EXISTENCE AND 
DISAPPEARANCE OF THE CONICAL 
SINGULARITY 


for the field at the origin, that is (j>'{r = 0 ) = 0 (stan¬ 
dard boundary condition), (j){r = 0 ) = (j)c (non-standard 
additional boundary condition). Reflecting the fact that 
the field equation of motion is of second order, these two 
boundary conditions completely fix the dynamics of the 
field even at infinity. This would make the system over¬ 
constrained, in general, and would reduce the freedom of 
solutions. 

The second possibility, which will be explored in the 
following, consists of imposing that, in a Taylor expan¬ 
sion of an around the origin, the constant term Q!h(^c, 0 ) 
identically vanishes for any (pc- In general this condition 
cannot be imposed by any symmetry, so that we need 
to tune the choice of the function an- Moreover, it is 
not clear whether such a functional form of an is valid 
in the presence of quantum corrections and whether such 
corrections spoil (or not) the classical picture by making 
the conical singularity reappear. The analysis of this im¬ 
portant issue is beyond the scope of our paper, but we 
would like to discuss it in a future project. 

To explore the second possibility further, we consider a 
canonical scalar field p without the potential, i.e., A 2 = 
—A/2 and C 3 = 0. To go beyond the Horndeski domain, 
we need to choose functions A 4 and B 4 which violate the 
condition A 4 = 2 XB 4 X — B 4 . One way is to extend 
the scalar-curvature coupling appearing in Brans-Dicke 
theories [1^ or dilation gravity [3l| to the forms A 4 = 
—Mp[Fi ((/))/2 and B 4 = Mp^F 2 {(j))/ 2 , where Fi{(j)) and 
^ 2 ( 0 ) are functions with respect to p. Another way is 
to take into account the X dependence in the functions 
A 4 and i? 4 , as it happens for covariant Galileons [s^ and 
extended Galileons [^ . 

To accommodate these two cases, we focus on the mod¬ 
els characterized by 


The results in Sec. Hill and m have been derived un¬ 
der the conditions that (i) (/'(r) = 0 for r > 0 and (ii) 
an = constant. Since the boundary condition <^'(0) = 0 
is generic for the regularity of solutions [H, [ 2 ^ , this is 
consistent with the assumption (i) around r = 0. The 
remaining possibility for eliminating the conical singu¬ 
larity is to construct models in which the parameter an 
approaches 0 as r —>• 0 . 

Let us consider this second possibility in detail. In gen¬ 
eral, around the center of a compact object with spheri¬ 
cal symmetry and staticity, we should expect that A —>■ 0 
for regularity. Therefore, the parameter an may be ex¬ 
pressed in the form 

an = ani(pc,0) + OiX). (5.1) 

We need to impose q;h(/'c,0) = 0 for removing the con¬ 
ical singularity. To fulfill this condition, there are two 
possibilities: 1 ) to solve an{$c, 0 ) = 0 for pc if a solution 
(or even more than one) exists; 2 ) to impose that the first 
term in the Taylor expansion (i.e., the constant term) in 
A identically vanishes. 

The first possibility cannot be viable in general, as this 
corresponds to setting two different boundary conditions 


A 2 — ~2^ ’ C3 — 0 , 

A4 = -iM2iEi (</) + /! (A), 

B 4 = ^M^,F 2 {P) + f 2 iX), (5.2) 

where /i(A) and f 2 (X) are functions of A. From 
Eq. (12.7p the function A 3 is given by 

A 3 = M^iVxF 2 AP) ■ (5.3) 

The parameters an and at read 


an 


at 


1 


Ki 

-fiFt 


F 2 ) - (/l + /2 - 2 Xf2,x) 


,(5.4) 


an - 


2Xf2,x 

A 4 


(5.5) 


The difference between Fi{p) and F 2 {p) gives rise to a 
non-zero contribution to an- The term fi + f 2 — ‘2Xf 2 ,x 
in Eq. 63) does not generally vanish for theories be¬ 
yond the Horndeski domain. The covariant Galileon cor¬ 
responds to /i(A) = a 4 A^ and / 2 (A) = 64 A^ with 
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the particular relation 04 = 864 , in which case 

h + f2- 2Xh,x = (a4 - = 0. 

In Eq. (15.211 the constant terms which can exist in 
the functions fi{X) (where i = 1,2) are assumed to 
vanish. In other words, after the Taylor expansions 
of the functions fi{X) around X = Q, i.e., fi{X) = 
/i(0) + f-{0)X + 0{X'^), the constants /i(0) have been 
absorbed to the functions Fi{(j)). As we will see below, the 
theories with Fi{cj)) ^ F 2 {cj)) are plagued by the existence 
of the conical singularity. To avoid this problem, we need 
to choose the two constants appearing in —A 4 and B 4 are 
exactly the same as each other. This corresponds to the 
aforementioned tuning of functions required to eliminate 
the conical singularity at r = 0 . 

A. Theories with 7 ^ F2{4>) 


where the second term in C 4 survives only for m = 1 . 

Since the terms C 3 and C 4 in Eq. (15.81) approach con¬ 
stants as r —>■ 0 , we can absorb these terms into the cos¬ 
mological constant A appearing in A 2 . Then, in the limit 
that r —0, Eqs. (12.111) and (12.121) reduce to the same 
forms as those studied in Sec. lYl with non-vanishing 
constants an = ctt = F 2 {(j)c)/Fi{(j)c) — 1. Hence the con¬ 
ical singularity is present at r = 0 for the theories with 
Fi{cj)) ^ F 2 {cj)). 

We note that the non-analytic function (j)'{r) = 

(0 < p < 1 ) satisfies the boundary condition (/>'( 0 ) = 
0 , but the second field derivative diverges for r —>• 0 . 
This gives rise to the divergent terms in C 3 and 

C 4 , whose functional dependence is different from that of 
the term ^A^atjr^ in Eq. (12.111) . This means that such 
non-analytic field derivatives do not eliminate the conical 
singularity. 


Let us first consider the theories with Fi{(j)) ^ F 2 { 4 >) 
and /i(A) = / 2 (A) = 0. In this case we have an = at = 
F 2 {(j))/Fi{cj)) — 1. For the boundary condition (j)'{0) = 0, 
the field (j) approaches a constant value (j)c B.sr ^ 0. Then 
the parameters an and at behave as constants around 
the origin of a compact object, so the conical singularity 
cannot be avoided at r = 0 . 

Next, we proceed to the theories with Fi{(j)) ^ F 2 {(j)) 
and non-vanishing functions of fi {X ) and /2 (A). An ex¬ 
ample of functions /i(A) and / 2 (A), which encompasses 
both covariant and extended Galileons, is given by 


B. Theories with Fi{(j>) = F2(4>) 


We proceed to the theories satisfying Fi{(j)) = F 2 {(j)) 
with /i(A) and / 2 (A) given by Eq. (15.61) . In this case 
we have 


a4A'" -h (1 - 2 n) 64 A" 

Mp2 Ei(0)/2 - 04A™ 

a 4 A™ -h 64 A" 

Mp2 Ei(0)/2 - a4A- ■ 


(5.9) 

(5.10) 


/i(A) = a 4 A’", / 2 (A) = 64 A" , (5.6) 

where m, n are positive integers (to, n > 1 ), and 04 , 64 are 
constants. Since / 1+/2 — 2 A/ 2 , 7 f = a 4 A'"-|-(l — 2 n) 64 A” 
and -2Xf2,xlA4 = 4 & 4 nA"/[A/ 2 jFi((/>) - 2 a 4 A'"] in 
Eqs. (lOl) and (OD . the parameters an and at again re¬ 
duce to an = at = F 2 { 4 >c)/Fi{4>c) — 1 = constant around 
the origin for the boundary condition (j){r —> 0) = (j)c- 
This constant behavior of an and at corresponds to the 
case studied in Sec. IlYl 

We caution, however, that the result in Sec. m has 
been derived by assuming that all the terms Ci defined 
in Eq. (12.151) vanish. Let us consider the effect of these 
terms around the center of a compact object. To satisfy 
the regular boundary condition (/)^( 0 ) = 0 , we expand the 
field derivative of the form 

='^KprP , (5.7) 

p=i 

where Kp is a constant and p is a positive integer. We 
substitute Eq. (15.71) into Eq. (12.151) by using the functions 
given by Eqs. (15.21) and (15.3|) with Eq. (15.61) . under the 
condition that $ approaches a constant as r —>■ 0. It then 
follows that 

Cl —^ 0, C 2 —^ 0, 

C 3 ^ KiMp 2 , 

C 4 ^ -2KiM2je“2'*’Fi,0 -h 8 k? 046 -^“^ , (5.8) 


For the boundary condition 4’'{0) = 0, both an and at 
approach 0 in the limit that r —> 0 . 

For the field profile (15.71) around r = 0, the quantities 
Ci behave in the same way as Eq. (15.81) with the relation 
Fi = F 2 . Since C 3 and C 4 approach constants as r —0, 
the solutions to Eq. f2.11|) and (I2.12|) around the center 
of a compact body are of the same forms as Eqs. (14.21) 
and (lOl) . respectively, with an = 0 and A subject to 
modifications arising from constant C 3 and C 4 . Since the 
term an vanishes for r —>■ 0 , we can avoid the conical 
singularity at the origin. 

In the above argument the equation of (j) has not been 
solved explicitly, but we employed the Taylor expansion 
of the form (EZD with the boundary condition ^'( 0 ) = 0 . 
In Sec. m we shall derive solutions to the scalar-field 
equation of motion in concrete models under the approx¬ 
imation of weak gravity and show the existence of the 
field profile of the form E3) around the origin. Pro¬ 
vided that the conical singularity is absent, the solutions 
derived under the weak-gravity approximation should not 
cause any divergent behavior for curvature quantities. 

VI. MODELS WITHOUT THE CONICAL 

SINGULARITY (Fi = F 2 ) AND VAINSHTEIN 
MECHANISM 

In this section, we solve the equations of motion for 
the field (j) and gravitational potentials 4>, 'k on the 
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weak gravitational background characterized by |$| <C 
1 , 1 'I'I <C 1 in concrete models without the conical singu¬ 
larity (Fi = F 2 ). 

We consider the situation in which the dominant con¬ 
tributions to Eqs. (I2.11|) and (12.121) are of the orders of 
and to recover the general relativistic 

behavior in the solar system (see Refs, [l^, [ 2 ^ for the 
similar approximation). From Eq. (12.141) the pressure Pm 
of non-relativistic matter is of the same order as 
so Pm is second-order in 'k. We deal with other terms 
(including an and at) as next-order corrections to the 
leading-order terms. 

Eliminating the terms $' and 'k' in Eq. (12.131) by us¬ 
ing Eqs. (j2.11l) and ()2.12|) . we obtain the equation for 
□'k coupled to □(/), where □ = cP/dr^ + {2/r){d/dr). 
Differentiating Eq. (I2.12p with respect to r and substi¬ 
tuting it into Eq. (I2.14L we can derive the equation of 
motion containing the terms □'k and □</>. Eliminating 
the term □'k by using these two equations, the resulting 
scalar-field equation of motion reads 

□ ~ piPm + P2 , (6.1) 


where 


d-i = 


M2 = 


and 


4’'r{A3^x +/3) 


4^4/3 

A 2 ,(j> 


( 6 . 2 ) 


— XA2^(I,x ) + (^ 3,0 — ‘^XA3^^x 


1 

fdr 

+A(j)'XA2^xx)r + ‘24''{^3 ,x + '^XA^^xx) — 

ji'XAi^xx + rai — A(j)'a2 


-\- 2 XA^^(f,x — 


, (6.3) 


where 9 , 04,64 are constants and n (> 1 ) is a positive 
integer. Then, the quantities pi and ^2 appearing in 
Eq. (16.11) are given, respectively, by 


Ml 


M2 


{qMpiF - jdpy 
2 p{Ml,F - 2040'"") ’ 

4(04 — 64 )n( 2 n — 1) 0'"" ^ 
0 P 


where 


0 



1 -I- 4(04 — b4)n{2n 


1 ) 


0/2(n-l) 

j,2 


( 6 . 10 ) 

( 6 . 11 ) 


( 6 . 12 ) 


A. Field profile and gravitational potentials around 
the origin 

We derive the solution to Eq. (EH) around the center 
of a compact body whose density approaches a constant 
Pm as r —>• 0. We search for regular solutions described 
by Eq. (15.71) around r = 0 for non-zero values of q. 

When n > 1 the terms 00' and 2040 ' " approach 0 
in the limit that r —>■ 0 , so the quantity pi reduces to 
Pi —>■ —qr/{2Mp\(3). Integrating Eq. (16.ip with respect 
to r, we obtain the following implicit solution 

r"0'-b 4 n(a 4 - 64 ) 0 '"""^ = , (6.13) 

oMpi 

where the integration constant is set to 0 to satisfy the 
boundary condition 0'(O) = 0. From Eq. (16.131) it is clear 
that there is a dominant solution of the form 


= i^2,x + ‘2XA2^xx)r + 2(A3 ^x + “^XA^^xx) 
'iXAi^xx , 2a2 


(6.4) 

(6.5) 

4,x ■ (6.6) 

Here and in the following, the kinetic term X should be 
understood as 0'". Under the above scheme of approxi¬ 
mation, Eqs. (12.111) and (I2.12p reduce, respectively, to 


r r 

ai = 2 XB 4 ^^x — -64,0 — A^^ 
a 2 = 2 XB 4 XX + B^^x — A 


/ Pm' 


4 A 4 

40' 


r 

4 A 4 


(^ 4,0 + 20"A4_jif 


A 2 — + 20"H3_x) 

1 




(6.7) 


^f' ~- 


r 

4 A 4 


(4,-2,-4,,) + ^ 


0'"A4. 


X 


A^r 


an 

2 r 


( 6 . 8 ) 


For concreteness, we study the model described by the 
Lagrangian (15.21) with the functions 

Fi(0) = F2(0) = F(0) = , 

h{X) = 04 ^" , / 2 (A) = 64 A" , (6.9) 


0 '(r) = cr, 


(6.14) 


where c is a constant. The coefficient c is different de¬ 
pending on the values of n, i.e.. 


c -b 8(04 — 64)0^ = 

qpm 


qpn 


C = 


3M 


pi 


3Mpi 
(for n > 2). 


(for n = 2), (6.15) 
(6.16) 


When n = 1, the solution to the field equation can 
be written in terms of an error function. Expanding it 
around r = 0 and using the boundary condition 0'(O) = 
0 , the dominant solution is given by 


0 '(r) 


_ qpm _ ^3 

12Mpi(a4 — 64) 


(6.17) 


Since both Eqs. (16.141) and (I6.17P are of the form (15.7L 
the discussion in Sec. IV Bl to show the absence of the 
conical singularity is valid. 

Integrating the solution 0'(r) = cr^ (where p = 1 for 
n > 2 and p = 3 for n = I), it follows that 0(r) = 
0c-bcrP+^/(p-b I) with (j)c an integration constant. Then, 
in the limit that r —> 0, the function F{(j)) converges to a 


















































constant value Fc = F((j)c). In this case the parameters 
q;h and at are in proportion to so they vanish at 

the origin. Employing these results and picking up the 
dominant contributions to Eqs. (IgTl) and (16.81) around 
r = 0 , they are integrated to give 

^ Pmr^ qcrP~^^ (1 + 

^ QMl,F,~ Mpi + 4(2p + 3)Mp2E, 

[(4pn+ 1)04 + b 4 ](crP)^" c$ 

( 2 pn+l)M 2 jFc r ’ 




Pmr'^ 

uuy, 

, [1 + (2 


qcrP+^ ^ (p + 2 + 

(1 + p)Mpi ^ 4(p + l)(2p + 3)M2 iE, 


n — l)p](a 4 — 64 )(cr-P)^" 
p{2pn + l)M^^Fc 


—+c^,(6.19) 
r 


where 0 $ and c^s, are integration constants. We need to 
choose c$ = 0 for the regularity of $ at r = 0. The 
constant , which can be derived by matching solutions 
at the radius of a compact body, is finite as in the case 
of the interior Schwarzschild solution. The contributions 
coming from at and an give rise to terms proportional 
to in Eqs. (I6.18|) and (I6.19L which vanish as r —>■ 0. 

Evaluating the Ricci scalar R by using the solutions 
(16.181) and (16.191) with 0 $ = 0, it follows that R ap¬ 
proaches a finite constant as r —>■ 0. Thus, the model 
(16.91) is free from the problem of the conical singularity 
as a consequence of vanishing parameters an and at at 
r = 0 . 


B. Vainshtein mechanism outside the compact 
body 

Since there is no conical singularity for the model 
(16.91) at the center of a compact object, we proceed to 
the derivation of solutions outside the body. In covari¬ 
ant Galileons [s^ with the dilatonic coupling F{<f>) = 
g- 2 <j<>/Mpi^ it is known that the presence of the terms 
/i(X) = 04 and f 2 iX) = 64 ^^ with 04 = 864 leads 
to the recovery of GR inside the so-called Vainshtein ra¬ 
dius rv, even for |g| of the order of 1 [H, [H, [s^- Now, 
we are going to discuss the Vainshtein mechanism for 
the model (16.91) without imposing the Horndeski condi¬ 
tion 04 -I- (1 — 271)64 = 0. Eor simplicity, we focus on the 
theory with n = 2 and 04 ^ 864 . 

The Vainshtein radius is characterized by the distance 
at which the field derivative in Eq. (16.121) becomes com¬ 
parable to the first term in Eq. (16.121) . i.e., 

~ 24104 - , ( 6 . 20 ) 

where ry can be explicitly known by solving the field 
equation (ICTl) . For the distance r ry the field self¬ 
interaction term p ,2 is suppressed relative to the term 
PiPm, with pi ~ qjMpi and /3 ~ —r/2 [I^. Then, 
Eq. (16.11) is integrated to give 


where Vg is the Schwarzschild radius of the source defined 
by 

Pm{r)f'^df. (6.22) 

Substituting the solution (16.211) into Eq. (I6.20p , we obtain 
the Vainshtein radius 

W = M^(24|o4-64|)-'/® , 

(6.28) 

where M is a constant having a dimension of mass. 

For Tg <C r <C cy the field-derivative term p 2 corre¬ 
sponds to the dominant contribution to Eq. dSH), such 
that piPm ^ P '2 — 2(/)'/r. Since the integrated solution 
in this regime is (j)'{ p) = constant, matching this solution 
with Eq. ()6.21() at r- = ry gives 



(j)'{r) = ^ for rg <C r < ry , (6.24) 

which is of the same form as that derived in Ref. [l^ in 
Horndeski theories. This reflects the fact that, provided 
04 7 ^: 64 , the factor p 2 is independent of the values 04 , 64 
and n. Hence the Vainshtein mechanism can be at work 
outside the compact body for the model (16.91) satisfying 
the condition 04 7 ^ 64 . 

The integrated solution to Eq. ()6.24p is given by ^(r) = 
4>q + qMpiVgr/vy, where (j)Q is a constant. Since the sec¬ 
ond term in (j)(r) is much smaller than Mpi for jgl < 1 and 
Tg <Sir <^ry, the quantity F[4>) = is close to 

1 (i.e., the value of GR) for |())o| <C Mpi. In the following 
we shall focus on the case |())o| Mpi and employ the 
approximation ~ —(2/Mpj)(l -|- 2 a 4 V^/Mpj). In¬ 
tegrating Eqs. dfiJ]) and (16.81) after substitution of the 
solution (I6.24|) . we obtain 


(j) ~ 


2 r 


1 - 2 q^ 


ry 


9^(1+ 8 g 2 ) Tg 


6 


ry \ry 


- 2(04 -I- 64 )— 2-3 —-r 


(6.25) 




2 r 


1 - 2 q^ 


_ 9^(1 -b 2g^) rg / 


ry J 


— 8(04 — 64 ) 


Mlyq^rl 


rln — 
r* 


ry \ry J 

(6.26) 


where r* is an integration constant. We have dropped the 
contribution coming from the term 2 a 4 X^/M^i in 
as this gives rise to corrections much smaller than the last 
terms of Eqs. (16.251) and (j6.26() . For the distance r <C ry 
the contributions other than the first terms in the square 
brackets of Eqs. (16.251) and (16.261) are much smaller than 
1, under which $ and —4' are close to the value rgj(2r) 
of GR. 

To quantify the deviation from GR, we define the post- 
Newtonian parameter 


b'{r) = 


qMpiVg 


$ 


7 = - 


4- ■ 


for r ^ ry 


( 6 . 21 ) 


(6.27) 
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The solar-system experiments placed the bound [13 

|7- 1| < 2.3 X 10“^ (6.28) 

On using the solutions (I6.25|) and (j6.26() . it follows that 


~ 1 


q^(l 4g^) 

2 ry 


2M^^q^r^r 


'V 


04 -|- 64 — 4(04 — 64 ) In ■ 


(6.29) 


The Vainshtein radius ry depends on the mass scale 
M defined in Eq. (16.231) . If the model (16.91) is responsible 
for the late-time cosmic acceleration, the parameter 04 
is related to the today’s Hubble radius rn ~ 10^® cm as 
l® 4 | dll' Assuming that 104 — 64 ! is of the 

order of I 04 I, we have ry ^ (j^jr'g from Eq. (I6.23|) . 

Since Vg ~ 10^ cm for the Sun, we can estimate ry ^ 10^° 
cm for jgj ~ 1. Inside the solar system (r < 10^^ cm), 
the second term on the rhs of Eq. (16.291) is smaller than 
10“®®. For r* between Vg and ry the term Inr/r* is at 
most of the order of 10 , so the last term on the rhs of 
Eq. (I6.29P does not exceed the order of 10“®®. Hence 
the experimental bound (I6.28|) is well satisfied inside the 
solar system. 

On using the relation |a 4 — 364 1 ^ ^’^cl the 

solution (I6.24p . the parameter an inside the Vainshtein 
radius can be estimated as 


jaHj 




4 4 


(6.30) 


which is smaller than 10“^® for jgj < 1. Thus, the devi¬ 
ation from Horndeski theories is significantly suppressed 
for r <C ry. If the Vainshtein mechanism is not at work 
and the solution is given by Eq. (I6.2ip outside a compact 
object, we have that janj ~ and hence janj ex¬ 

ceeds the order of 1 for r < y^jglrnr^. For the Sun and 
|(jf| ~ 1 this condition translates to r < 10 ®® cm, so janj 
is much larger than 1 in the solar system. The above ar¬ 
guments show how the Vainshtein mechanism is efficient 
to suppress both the propagation of the fifth force and 
the deviation from Horndeski theories. 

We recall that, for n = 2, the solution to Eq. (ICT) 
around r = 0 is given by (ji^r) = cr, where c obeys the 
relation (16.151) . For the mass scale M of the order of 
M ~ |a 4 |“®A ~ term 8(04 — 64 ) 0 ® is the 

dominant contribution to the Ihs of Eq. (16.151) . so the 
solution around the origin is given by 


(j)'{r) 


(qPmM^V'^ 

V Mpi ) 


(6.31) 


where we have assumed g > 0. In fact, the ratio ^ = 
c/[ 8 |a 4 — 64 IC®] is of the order of ^ ^ [Po/( 9 Pm)]^'^^ 1 

for g 1 , where po ~ 10 “^® g/cm® is the cosmological 
density and pm 1 g/cm® is the mean density of Sun. 


The validity of the solution ()6.31l) is ensured around 
r = 0 , but we can extrapolate it to the surface of 
the compact body (radius tq) to estimate the order of 
4>'{rp). On using the relations ~ Pm'TQ/M'^^ and 
M = (gMpirg)®/®/ry the extrapolation of Eq. (16.311) 
gives (l)'{rQ) ~ qM^irg/ry, which is of the same order 
as Eq. ()6.24l) . Thus, the two solutions (16.241) and (16.311) 
smoothly match each other around r = rp. There are cor¬ 
rections to the solution (16.311) , but they do not change the 
order of 4>'{r) inside the body. This situation is analogous 
to what was found in Ref. [Q in Horndeski theories. 

The smallness of an in our model comes from the fact 
that an depends on X alone (an oc V^). Under the 
operation of the Vainshtein mechanism, the suppression 
of the field derivative leads to the small value of an- 
This is in contrast with the models where an depends on 
(j) (i.e., ^ F 2 {(j))). In the latter case, the models 

suffer from not only the conical singularity problem, but 
also the breakingof the Vainshtein mechanism occurs as 
shown in Refs. 2ll-[^. The fact that Ei((/)) and F 2 { 4 >) 
are equivalent to each other is crucial for the success of 
the Vainshtein mechanism. 


VII. CONCLUSIONS 

In GLPV theories where the deviation from Horndeski 
theories is weighed by the parameter an, we have shown 
that the conical singularity arises at the origin of a spheri¬ 
cally symmetric body for nonzero constant an around the 
origin. For both vacuum and interior Schwarzschild so¬ 
lutions satisfying the boundary condition (j>'{r = 0 ) = 0 , 
the Ricci scalar diverges as R = —2an/r'^ around r = 0. 
In such cases, the spherically symmetric static body does 
not exist due to the singularity problem at its center. 
This divergence originates from the nonvanishing contri¬ 
bution Q!h to the gravitational potentials, which does not 
allow for the recovery of the Minkowski space-time. 

For the theories in which an vanishes as r —>■ 0, it is 
possible to avoid the appearance of the conical singular¬ 
ity. This requires the condition that the functions Fi {(f>) 
and F 2 { 4 >) appearing in Eq. (15.21) are equivalent to each 
other. The functions /i(A) and / 2 (A) need to be chosen 
such that they do not involve arbitrary constants which 
give rise to the difference between —A 4 and B 4 . Violation 
of the condition Fi ((/>) = F 2 (</) means that the geometric 
structure of the four-dimensional Ricci scalar R is mod¬ 
ified, which causes a geodesically incomplete space-time 
at r = 0 . 

The model described by the functions (16.91) , which cor¬ 
responds to the extension of covariant Galileons with a 
dilatonic coupling, is free from the problem of the con¬ 
ical singularity. This model is designed to have the de¬ 
pendence of an proportional to A" around r = 0. We 
derived the field profile as well as the gravitational po¬ 
tentials around the center of a compact object under the 
approximation of weak gravity and showed that the Ricci 
scalar remains finite as a consequence of the vanishing an 
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at the origin. 

For the model (|6.9I) with n = 2 we also found that 
the Vainshtein mechanism is at work outside the body 
to suppress the fifth force and the parameter oh in the 
solar system. The regular solution of the field derivative 
4 >'{r) around the origin can be extrapolated to match the 
exterior solution around the surface of the body. 

There are several issues we have not addressed in this 
paper. First, it will be of interest to study whether the 
similar properties for appearance and disappearance of 
conical singularities also hold for the theories involving 
the Lagrangian L 5 related to the Einstein tensor. Sec¬ 
ond, we showed that the functional form of Lagrangian 
must be properly chosen for eliminating the conical sin¬ 
gularity, but the absence of a symmetry may give rise to 
quantum corrections which can modify the Lagrangian 
structure. If quantum corrections preserve the structures 
of four-dimensional curvature quantities like the Ricci 


scalar R, we expect that the modification is less harmful. 
If any difference of a constant between the functions —A 4 
and B 4 arises by quantum corrections, this would lead to 
reappearance of the conical singularity. We leave these 
interesting issues for future works. 
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